Mott transition at large orbital degeneracy: dynamical mean-field theory 
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We study analytically the Mott transition of the iV-orbital Hubbard model using dynamical 
mean-field theory and a low-energy projection onto an effective Kondo model. It is demonstrated 
that the critical interaction at which the insulator appears (U c i) and the one at which the metal 
becomes unstable (U C 2) have different dependence on the number of orbitals as the latter becomes 
large: U c i oc \/~N while U C 2 oc N . An exact analytical determination of the critical coupling U C 2/N 
is obtained in the large-iV limit. The metallic solution close to this critical coupling has many 
similarities at low-energy with the results of slave boson approximations, to which a comparison 
is made. We also discuss how the critical temperature associated with the Mott critical endpoint 
depends on the number of orbitals. 



I. INTRODUCTION 

The Mott transition is one of the central problems in 
the field of strongly correlated electronic systems. It 
is directly relevant to many compounds, such as V2O3, 
NiS2- £! ;Se 2 ; [Q, fullerenes [0, and two-dimensional organ- 
ics of the k-BEDT family ||. Viewed from a more gen- 
eral perspective, the Mott transition confronts us with 
the need to describe electrons in a solid in an intermedi- 
ate coupling regime, where perturbative methods around 
either the localised or the itinerant limit are of very lim- 
ited use. The Mott phenomenon is also a key issue for 
the development of electronic structure calculations of 
strongly correlated materials, a subject of current intense 
investigation 0. 

A breakthrough in the understanding of this problem 
resulted from the development of dynamical mean-field 
theory (DMFT) §, § |, [|. Following the observa- 
tion that a Mott transition as a function of the inter- 
action strength U takes place in a (fully frustrated) one- 
band Hubbard model at half- filling [|| within DMFT, 
a rather surprising scenario was unravelled JTcj ] . In the 
parameter space defined by temperature and interaction 
strength, two kinds of solutions of the DMFT equations 
were found, corresponding to an insulator and a metal, 
with a region U C \{T) < U < U C 2{T) in which both types 
of solution coexist. The spectral function of the insulat- 
ing solution has only high-energy features corresponding 
to Hubbard bands, while the metallic solution also has a 
low-energy quasiparticle resonance. Coexistence results 
in a first-order phase transition at finite temperature. 
The first order line ends into a second-order critical end- 
point at T = Tmit, U = Umit at which the coexistence 
domain closes (U c1 (T M it) = U c2 (T M it) = U M it)- At 
T = 0, the transition is also second-order, because the 
metallic solution always has lower energy than the in- 



sulating one. A qualitative picture of the evolution of 
the spectral function as a function of control parameters 
was obtained, and this resulted in successful experimen- 
tal predictions jll] . 

Early work on these issues relied heavily on an accurate 
but approximate technique, the iterative perturbation 
theory || 10 1. Later on, exact results near U c n(T = 0) 
were obtained (lj, thanks to the development of a 
renormalization method exploiting the separation of en- 
ergy scales that holds true close to this critical coupling. 
This method performs a projection onto a low-energy ef- 
fective theory after elimination of the degrees of freedom 
corresponding to Hubbard bands. A clear picture of the 
structure of the DMFT equations near the Mott tran- 
sition and of the different critical points was finally ob- 
tained through the introduction of a Landau functional 
approach [[L4|. These established results for the one- 
band Hubbard model together with numerical methods 
for solving DMFT equations, such as Quantum Monte 
Carlo [ll| and Numerical Renormalization Group [|l6| , 
provide now a consistent picture of the Mott transition 
within DMFT. 

In this paper, we investigate the Mott transition in the 
limit of large orbital degeneracy, within DMFT. We es- 
tablish analytically that the critical couplings U C \{T = 0) 
and U C 2(T = 0) are not only different, but that they 
have very different scalings as a function of the number 
N of orbitals. Indeed, U c i increases as y/W, while U C 2 
increases as N, so that the coexistence region widens. 
Close to U C 2, the separation of scales which is at the 
heart of the projective method becomes asymptotically 
exact for large N. This allows us to determine exactly 
the value of the critical coupling U c2 /N at large N, and 
also to demonstrate that slave-boson like approximations 
become asymptotically valid at low energy in this limit. 
Finally, we discuss how the finite temperature metal- 
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insulator transition depends on the number of orbitals. 

Our findings on the critical couplings are in agree- 
ment with early QMC results for the two-orbital model 
(N = 4) [M , and with recent results for higher values of 
N p^j . A widening of the coexistence window [U& , U C 2\ 
as N increases is clearly seen in these simulations. Also, 
for fixed N , the critical coupling required to enter the 
Mott state is largest |Q at the particle-hole symmetric 
filling n = 1/2. The present work also puts the results 
of Refs.|^9|, in a new perspective: there, a VN scal- 
ing of the critical coupling was proposed. We find that 
this indeed applies to the coupling where the insulator 
becomes unstable (£/ c i), while the true T = Mott tran- 
sition (at which the quasiparticle residue vanishes) takes 
place at U C 2 oc N. Indeed, in Ref.|l?J, the y/N scal- 
ing was rationalized on the basis of a stability argument 
for the insulator. Finally, the rather high temperatures 
considered in |20|j explain why only the y/N dependence 
was reported there: distinguishing U C 2 from U c i requires 
significantly lower temperatures fl22| . 

The limit of large orbital degeneracy is relevant to the 
physics of systems with partially filled d- or f- shells, 
as well as to fullcrcncs. Direct numerical approaches 
become prohibitively difficult as the number of orbital 
increases (QMC methods scale as a power law of N, 
while exact diagonalizations |l7|] scale exponentially). It 
is therefore important for future research to develop ap- 
proximate but accurate impurity solvers which can han- 
dle many orbitals. The controlled results that we estab- 
lish in this paper can be used as tests of these approxi- 
mation methods. 



II. MULTI-ORBITAL HUBBARD MODEL 

We consider a generalized Hubbard model involving N 
species of electrons, with Hamiltonian: 
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where i,j are sites indices, a the orbital index, Jl is the 
chemical potential and n the average density of particles 
per species : 



^£<4* 



(2) 



Introducing n in the Hamiltonian is just a convention for 
the chemical potential Jl. In particular, it is convenient 
at half-filling where the particle-hole symmetry implies 
n = h and Jl — 0. For a single site (atom), the spectrum 
consists of N+\ levels, with energies ^{Q—nN) 2 depend- 
ing only on the total charge on the orbital: Q = 0, • • • ,N 
and with degeneracies (g) . 



The usual single-orbital Hubbard model corresponds to 
N = 2 (cr =t,l)- The Hamiltonian considered here has 
a full SU(A r ) symmetry which includes both spin and or- 
bital degrees of freedom. Starting from a more realistic 
model which assumes an interaction matrix U mm i be- 
tween opposite spins and U m m' — Jmm' between parallel 
spins, it corresponds to the limit of isotropic U m m' = U 
and vanishing Hund's coupling J mm i = 0. 

When U is large enough, we expect a Mott insulating 
state to exist at fillings n — Q/N, corresponding to an in- 
teger occupancy of each site on average Q = 1, - • • ,N— 1. 
In this paper, we investigate analytically the nature of 
these Mott transitions within DMFT. We consider only 
phases with no magnetic ordering and study the transi- 
tion between a paramagnetic metal and a paramagnetic 
Mott insulator. In Section p| , we extract the large-iV 
behaviour of U c \ and U C 2 using the low-energy projective 
technique analysis of the DMFT equations introduced in 
[IjJI for N = 2 and extended in [M : we show that the 
equations for the t/ c 's derived by this method are greatly 
simplified for N — > oo in the sense that an atomic limit 
becomes exact. In section |lV| we find quantitative agree- 
ment between these results and a multiorbital slave boson 
method. Finally, in Sec.|v| we consider the Mott transi- 
tion at finite temperature. 



III. LARGE- iV BEHAVIOUR OF U cl AND U c2 

A. DMFT and the low-energy projective method 

Let us recall briefly the DMFT equations and their 
low-energy projective analysis |l2|, O. DMFT maps the 
lattice Hamiltonian above onto a multi-orbital Anderson 
impurity model with the same local interaction term than 
in (Q) and a hybridization function A(iu) n ). The local 
Green's function reads: Gd(?w„) _1 = iiv n + Jl — A(iu n ) — 
In this expression, Jl is the chemical potential 
(shifted by the Hartree contribution) and E^ia-Vi) is a 
local self-energy. A self-consistency requirement is im- 
posed, which identifies Gd(iw n ) with the on-site Green's 
function of the lattice model with the same self-energy 
Sd(iw n ). This reads: 



G d (iuj n ) = de 



D(e) 



S d (iw„) 



(3) 



In this expression, D(e) is the free electron density of 
state (d.o.s.) corresponding to the Fourier transform of 
the hopping matrix elements ty . In the particular case 
of a semi-circular d.o.s. with half-width D = 2i, equa- 
tion (j^) simplifies to: 



A(iu> n ) = t 2 G d (iu) n ) 



(4) 



Solving the impurity model subject to (Q) determines 
both the hybridization function and the local Green's 
function in a self-consistent manner. In order to give 
an explicit Hamiltonian form to the Anderson impurity 
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model, the function A can be represented with an aux- 
iliary bath of conduction electrons Cfc CT , with effective 
single-particle energies e k and effective hybridizations Vk 
to the local orbital such that : 



AM = Y 



\v k \ 2 

iui - e k 



H = Y ^ C \,a C K<y + Y V *( C k,<T d <> + d Uk,a) 



u 



X)(4^-») 



(6) 



Within DMFT, the Mott transition at U c2 and T = 
is associated with a separation of energy scales. For U 
slightly smaller than [/ c2 , the d-electron spectral function 
has a sharp quasiparticle peak, well separated from the 
high-energy Hubbard bands. Naturally, a small amount 
of spectral weight does connect these two features in the 
metal. Slightly above £/ c2 , the Mott gap A g separating 
the Hubbard bands of the insulating solution is finite. 
While the weight Z of the quasiparticle peak vanishes 
at U C 2, the Mott gap A g disappears only at the smaller 
critical coupling U c \. Because of the self-consistency con- 
dition (^), the hybridization function of the effective con- 
duction bath (and hence the e^'s) also displays this sep- 
aration of scales. 

This remark led the authors of Ref. jl^] to use a pro- 
jective method in order to study the critical behaviour at 
U C 2- Within this approach, the effective conduction elec- 
tron bath is separated into a high energy sector (k £ H ) 
associated with the upper (H + ) and lower (H~ ) Hubbard 
bands, and a low-energy sector (fc S L) associated with 
the quasiparticle resonance. The impurity model Hamil- 
tonian is thus written as (for simplicity we take Vl and 
Vh independent of k): 



H — i?L + #H + #HL 

#l = Y efec L c ^ 



(7a) 
(7b) 



fcSL.er 



Hr= Y [ e kc\,a c K<y + ^ c k,a d <r + h.c 



fceH±,o- 



u 

~2 



Y d * d ° - n 



- v Y d * d ° 



hul = Y [ v ^ c L d ° + h - c ] 



(7c) 
(7d) 



fcSL.cr 



We note that the self-consistency condition (||) implies 
that J2keL^L A 2 ^ s proportional to the quasi-particle 
residue Z, so that it is clear that the low-energy hy- 
bridization vanishes at U C 2- 

In order to obtain an effective theory close to U C 2, one 
first diagonalizes the high-energy part H-r of the Hamil- 
tonian, and obtains an effective low-energy hamiltonian 
by expanding in Vl- From the presence of a finite gap in 



the high-energy bath (for k € H 1 * 1 ), Hn describes an 
impurity in a semiconducting bath. The position of the 
d-level (i.e. the chemical potential) fixes the number of 
electrons in the ground-state. For a charge Q — nN, this 
high-energy problem has a (q) degenerate ground-state, 
(5) spanned by eigenvectors (with eigen-energy E gs ). The 
detailed form of these states is complicated (they mix 
with higher impurity charge states through Vh), but they 
form a representation of the SU(N) "spin" group, of di- 
mension (q)- 

Using a generalized Schrieffer- Wolff (SW) transforma- 
tion to order in order to project onto the low-energy 
Hilbert space (defined by the selected charge sector of the 
impurity d, and the low-energy sector for the conduction 
bath k G L) , one obtains an effective Kondo Hamiltonian 

§§§: 



Y 

fee L,cr 



U ^ " lir 

kk' .jiv : aa 



(8) 



where X** v = is a Hubbard operator connecting 

two states in the ground state manifold. The Hamilto- 
nian contains potential scattering terms for a = a' and 
spin- flip (Kondo) terms for a ^ a' , with the spin opera- 
tor S^v = X^ v — S^X" . The matrix of (dimensionless) 
coupling constants reads O, |23| : 



T(T(T 

J p,v — 



u 



fi 



U 



-dl 



(9) 



We emphasize that the low-energy Hamiltonian has been 
derived to first order in V^: this is asymptotically jus- 
tified exactly at the critical point (and hence sufficient 
to determine the critical coupling), but deviations away 
from the critical point require to consider higher order 
terms in Vl . Finally, we note that the bandwidth of this 
Kondo problem is of order Z t, while the coupling is of or- 
der V^/U ~ Zt 2 /U. Hence, the dimensionless coupling 
constant of the effective low-energy Kondo problem is 
of order one: we have an intermediate coupling Kondo 
problem [ jtjj . 

The conditions that determine t/ cl and U C 2 have been 
derived within the projective method Jl2| , p3[ and within 
the more general Landau functional approach |l4|] (in 
which the stability of the insulating solution to differ- 
ent kinds of low-energy perturbations is studied using a 
Landau functional of the hybridization function A(iu> n )). 
The condition for U C 2 reads: 



( N 



(10) 



{Xp,vC a c f)) ( Jpv Jpfi 



p,p °vp 



As explained in an analogous equation is found 

for f/ c i , where off-diagonal terms in the average are taken 
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to zero. For finite N, this equation is complicated to 
solve, since the average on the right-hand side is a func- 
tion of J which we can only obtain by solving the low 
energy Kondo problem explicitly. However, we will show 
in the following that : 

1. Using the SU(N) symmetry, the tensor structure 
of J can be parametrized with only two numbers 
A and B, which are dimensionlcss functions of N, 
n, an d Va/U. This greatly simplifies ( |To| ) (Section 

Si - 



2. An explicit calculation of A and B can be per- 
formed in the atomic limit Vh = 0. This allows 
us to show that the two critical couplings have the 
following TV-dependence at large N: 



U cl = Vnu c1 
U c2 = NU c2 + 



(11a) 
(lib) 



In fact, the atomic limit evaluation of A and B 
becomes asymptotically exact when U is taken to 
be proportional to N. This is not the case when 
U oc yN, but we show that corrections to the 



atomic limit do not modify the scaling (11a) 



3. For U oc N, the effective Kondo problem can actu- 
ally be solved explicitly at large N, using standard 
techniques. This allows for an explicit computa- 
tion of the matrix element (yX^c^c^ involved in 
(|l0|). Since, from the above remark, the effective 
couplings entering the Kondo problem can be eval- 
uated in the atomic limit in that case, this allows us 



to determine the prefactor U c2 exactly (Sec. |III D| ) . 
Because the atomic limit is not asymptotically ex- 
act for U oc */N, we are not able to obtain an exact 
determination of U c ±. 



B. Simplication due to the SU(N) symmetry 



with Q = nN boxes), the tensor structure of J is then 
given by the Wigner-Eckart theorem : we have to find 
the occurence of R in (id © ad) ® R. We have a triv- 
ial term id ® R = R and a single non trivial one since 
ad® R — R ® ... , where . . . denotes other irreducible 
representations. Then we have : 



1 



;t 



iT H -E 



fj, ) = A^aa'Sfj 



(13) 



where C are the Clebsch-Gordan coefficients (R\ad (g> R). 

In order to avoid the computation of C, we can use the 
following argument. Let us consider the operator d a d\,, 
which has the same symmetry as Oi, so the previous 
symmetry argument gives 



d„d 



/' 



Ai&cjcj'&iiv 



B 2 C° 



(14) 



with the same C, but different coefficients A 2 ,B 2 . Since 
B 2 =/= 0, we can eliminate C between the two equations. 
We finally obtain the general form, valid to first order in 
V 2 for a SU(N) symmetric model: 



T>6. 



B 



A 4 



dld B 



(15) 



In this expression, A (N, n; Vk/U) and B (N, n; Vr/U) 
are dimensionless coefficients which depend only on 
the dimensionless ratio Vn/U associated with the high- 
energy problem (and of course a priori also on orbital 
number N and density n). 

Inserting ( |l5| ) into equation ( |l0| ) , we obtain the simpli- 
fied form of the conditions which determine the critical 
couplings: 



IT 2 

t 2 
TI 2 

t 2 



: [A + nB] 2 + B 2 n(l-n)(N + 1) 
: (A + n B) 2 + B 2 n(l - n)(N + 1) 



(16a) 



(16b) 



Computing the couplings (g) is a complicated task in 
general, since it requires a knowledge of the states \u), 
and thus a diagonalization of Hn. However, general con- 
siderations based on the SU(TV) symmetry of the model 
allow one to reduce the unknown Kondo couplings J™ 
to two dimensionless numbers only. Let us consider the 
operator: 



or =d 



i 



-dl 



(12) 



Since the Hamiltonian is SU(iV) symmetric, this oper- 
ator transforms under SU(-/V) with the representation 
F (8 F* = id © ad ; the tensor product of the funda- 
mental F (one box in the Young tableau language) and 
its complex conjugate F* (N — 1 boxes), which reduces 
to the sum of the identity id and the adjoint ad. If we 
denote by R the ground state representation (a column 



C. The atomic limit 

The atomic limit approximation is defined by Vh = 0. 
In this limit, the states \n) are just the impurity states 
|/i)o = \<Ti ■ ■ ■ <tq) (Q = 0, • • • , N), without mixing with 
the high-energy electron sector. In that case, A and B 
(from (|l5|)) can be computed explicitly ( we choose here 
Jl = 0, i.e. the chemical potential is bound to be at the 
center of a step in the Coulomb staircase). We find: 



.4 



-2,B 



(17) 



We note that these coefficients do not depend on N in 
the atomic limit. Moreover, it is clear that, in general 
(for an explicit proof, see Scc4IID) : 



(18) 
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When using these results into (16a-16b), we obtai n the 
TV-dependence of the critical couplings quoted in (11a- 



11b) 



This derivation relies on the use of the atomic limit 
Va = 0. An important question is therefore to check that 
the TV-dependence of the critical couplings is not modi- 
fied by the interaction Vh that dresses the atomic states 
We examine carefully this question in Appendix [A|, 
where we give an explicit proof that this is indeed the 
case. The problem is in fact quite different for U c i and 
U C 2- For the latter, U must be scaled as U — NU. In 
that case, we show in the appendix that the atomic limit 
calculation of A and B is in fact exact at large TV. This 
can be expected since the band gap diverges with TV, so 
that the Hubbard bands are well separated from the low 
energy degrees of freedom. However, one might worry 
that the growing number of orbitals (TV) could counter- 
balance this effect. This is actually not the case, and we 
show that the corrections to J (or A,B) coming from Vh 
are subdominant in large-TV, order by order in pertur- 
bation theory. Moreover, the low-energy effective Kondo 
problem becomes exactly solvable when U oc TV (i. e clos e 
to U C 2)'- this is the subject of the next section (Sec. 1IID| ), 
in which the exact value of the prefactor U C 2 is obtained. 

The problem is more involved for U c \, corresponding 
to the scaling U ex VTV. In that case, corrections in Vh 
beyond the atomic limit do produce corrections to A and 
B (hence to J). However, we are able to show (Sec. A 2 
of the appendix) that these corrections are of order one, 
so that the scaling U& — VNU c i does hold. 



D. Explicit large- TV solution of the effective low 
energy Kondo problem at [7 C 2 



In this section, we explicitly solve the low energy effec- 
tive Kondo problem in the large- TV limit using standard 
methods |2(|, in order to obtain the prefactor U C 2- The 
effective Hamiltonian (||) reduces at order (1/TV 2 ) to: 



fceL.er 



(19) 



NU , 



ft f _ ° aa ' 
J a' J a 2 



r Ck> 



This model becomes exactly solvable at large TV because 
the Kondo coupling scales as 1/TV (which is a result of 
U oc TV). Introducing an auxiliary boson field b conju- 
gate to Ylka- fa Cka > an d a Lagrange multiplier field A to 
enforce the constraint flf<? — Nn, we see that, as 
N — > oo, the field b undergoes a Bose condensation, and 
the saddle-point values of the 6, A fields are determined 



by the equations (at T = 0): 
G7 1 H=c + A-6 2 G° L H 



Uci 



<\uAmG j(uj) 

3 

dTG f (T)G° cL (T) = 



(20a) 
(20b) 

dcIm[G / ( W )G° L (c)] 
(20c) 



In this expression, G° L denotes the Green's function of 
the low-energy part of the effective bath of conduction 
electrons, in the absence of Kondo coupling, namely: 



1 



keh 



y2 AlM 



(21) 



First we express the self-consistency condition (|3|) 
or more precisely its low-energy counterpart. We re- 
late the low-energy part of the on-site Green's function, 
Gdh{ioJ n ), to the auxiliary-fermion one, Gf(iu> n ), by im- 
posing that the low-energy part of the interacting Green's 
function of the effective bath of conduction electrons, 
G c l , can be calculated either from the original Anderson 
model or from the low-energy projected Kondo model, 
with identical results. G c l is related to G° L through the 
conduction electron T-matrix: 



G c l — G" L + (G° L ) 2 T 



(22) 



The Anderson model T-matrix is V^Gdh, while it is b 2 G/ 
for the low-energy Kondo model. Therefore Gdh is di- 
rectly related to G/ by: 



VjjGa = b 2 G f =T 



(23) 



Using this into Eq.(20a), one finds the expression of Gdh 
in terms of the low-energy hybridization Al: 



1 



(iuj n + A) - A L (iuj n ) 



where we have introduced the notation: 

z- h - 



(24) 



(25) 



We can finally obtain Gdh in closed form by making use 
of (|24|) into the self-consistency condition (|J) (restricted 
to the low-energy sector). This yields: 



G d h{iu) 



de 
de 



D(e) 



AlM + G^M 
D{e) 
+ A) — e 



(26) 
(27) 



The low-energy part of the local spectral density reads: 



p dL (u>) = — ImGdL(w + i0 + ) = D 

IT 



ui + A 



(28) 



6 



By comparing this form of Gdh to its expression in terms 
of the local self-energy, we recognize that Z is the quasi- 
particle weight and A yields the zero-frequency limit of 



M-S rf (0) = ^ , Z = 



1 5Srf l ' 



(29) 



The value of X/Z can b e det ermined by using the explicit 
form of pd = Zpf into (20b): 



d^ D (^±± 

, Z \ Z 



(30) 



Changing variables, this leads to X/Z = fXo(n) where 
fio(n) is the non-interacting value of the chemical poten- 
tial, defined by: 



/Mo 
deD(e) 
-CO 



(31) 



We note that this is precisely the result expected from 
Luttinger's theorem (we have p, — £^(0) = po( n ) which 
insures that the Fermi surface volume is unchanged by 
interactions). Hence, the low-energy part of the spectral 
function (quasi-particle resonance) is simply given by: 



p dh (uj) = D +/x (n)J 



(32) 



Its shape is simply given by the non-interacting d.o.s., 
centered at u) = —Zp and with a width renormalized by 
Z. 



Second, we can now use ( P0q ) and Im[AiG/] 



^ImG/, to obtain 



Uc2 

4 



1 

Z* 



duo(uj + A)pdi(w) 



(33) 



D { — ^- ) (34) 



Z 



z 



Clearly, Z drops out of this equation, which has there- 
fore a unique solution, the critical coupling U C 2 at which 
Z vanishes. This is expected: the low-energy effective 
Hamiltonian is valid only exactly at the critical point. A 
determination of the critical behaviour of Z slightly away 
from the critical point would require considering higher- 
order terms in the effective low-energy Hamiltonian. We 
therefore obtain the large- TV value of the critical coupling 
(using the expression of X/Z found above): 



with: 



U C 2 = ^=4e Q (n) 



dee D(e) 

-oo 



(35) 



(36) 



This result is reminiscent of the critical coupling obtained 
in the Gutzwiller or slave-boson methods |E3| : a detailed 



comparison is made in Sec.tV. Taking for example a 



rectangular d.o.s. of half-width D, Eq.(pq) yields: 

Uc2 

N 



ADn{\ - n) 



(37) 



U C 2 is largest at half- filling (n = 1/2). This is because 
orbital fluctuations are maximum there, so that it is more 
difficult to destroy the metallic phase (see however jl8| ) . 
Obviously, for finite N, a Mott transition is found only for 
rational values of n — Q /N with integer Q (these values 
are dense on the interval [0, 1] as N becomes large). 

Anot her w ay to obtain U C 2 is to use the criterion of 
section [II C together with the present large- N limit in 



order to compute the matrix element ^2 aa i ^S' (T(T 'cj T ,c c 

This can be exactly evaluated at N = oo using the effec- 
tive action for the Kondo model (|20|): 

S cS = [ dTVte-AV(/]/ ff -n) (38) 

rP NT J 

+ / dT—^bi+Y.bfW+bciu 

JO 4l/ L 

[ dr [ Ar'Y,cl{r)Gl lj {T-T')- 1 c a {T') 
Jo Jo a 

Because this action is gaussian, it is straightforward to 
get the following average (here a ^ a'): 



dTG° cL (r)G f (r) 



We find therefore: 



'C',C C 



-N z b 



2,2 



drG f (r)G cL (r) 



(39) 

O(N) 
(40) 



which indeed scales as TV 2 . Inserting this into (16b), we 
have : 



U, 



c2 



At 



Vl 



Zt 



Nb / drG f (r)G Q cL (T) (41) 



which coincides with Eq. (20c) since yZ = b/V^. Hence, 
we have checked that both methods (expl icit large- N so- 
lution and projective criterion of Sec. [III Q ) yield the same 
result for the critical interaction U C 2- 



IV. COMPARISON TO A SLAVE-BOSON 
APPROACH 

In this section, we show that a multiorbital slave bo- 
son approximation reproduces the previous results close 
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to the critical coupling U c %. This method is an extension 
to many orbitals of the Gutzwiller approximation, for- 
mulated in terms of slave bosons by Kotliar and Ruck- 
enstein E^j and extended to many orbitals in |2q] . In- 
deed, with the scaling U = NU, the Mott gap is of order 
N and quasiparticles dominate the physics over most of 
the energy range. Furthermore, the expression (|2g| ) for 
the physical spectral function close to U C 2 contains pre- 
cisely the same ingredients than in slave boson methods, 
namely: a shift of the effective level which insures Lut- 
tinger's theorem, and a quasiparticle weight Z (vanish- 
ing at the transition) . The self-energy contains no further 
renormalizations (in particular the quasi-particle lifetime 
is infinite). This bare-bone picture of a strongly corre- 
lated Fermi liquid is precisely that emphasized by slave 
boson approaches. 

Following pqj , one introduces 2 N slave bosons 4' ai ... C j p 
associated with a given spin configuration <j\ ■ ■ ■ a v on 
the impurity orbital. The probability that the or- 
bital is in a given spin configuration is simply given by 
(<%■!.. .a- <?W...o>)- The physical electron operator is rep- 
resented as: 



d a z a fa 



(42) 



with: 



91 



E E 



W..^^!...^ 92 (43) 



P a^(cr u ...,a p ) 



9l 



92 



1_ E E C.^Al..^ 

P ae(<Ti,—,<T p ) 



i-E E 



^a 1 ...(j„4 > (J 1 ...a p 



-1/2 



-1/2 



P a<£(ai,...,a p ) 

This representation is supplemented by two constraints: 

/i/« = E E C^-* ( 44a ) 

P ae(o-i,...,CT p ) 

1 = E E ( 44b ) 

P (<Ti...Cr p ) 

At this stage, this is an exact representation of the 
Hilbert space, which can be used to rewrite the Hubbard 
Hamiltonian. In order to obtain an approximate solu- 
tion, we shall assume that all auxiliary bosons undergo 
a Bosc-condcnsation, i.e become static c-numbers, and 
that furthermore the corresponding expectation value de- 
pends only on the total charge p and not on the specific 
spin configuration a±, ■ ■ ■ ,cr p . We shall thus set, at each 
site: 



"ia 1 ...a p — Yp 



(45) 



We emphasize that this approximation does not corre- 
spond to a controlled saddle-point, even in the large- iV 



limit, since the number of auxiliary fields grows rapidly 
with N. 

Using this Ansatz, the Hubbard Hamiltonian can be 
rewritten in the form: 



H = 



with: 



N 

■EE**. 



JJ N 

|EE 

i p=0 



4> 2 p {p-Nn) 2 
(46) 



JV-l 

E 

p=0 



N - 1 
P 



(47) 



and the constraints: 

N 



E 

p=0 



E 

p=l 



N - 1 
p-1 



n = {ft fa) 



(48) 
(49) 



The free-energy corresponding to this Hamiltonian must 
then be minimized with respect to the N + 1 variational 
parameters p 's, subject to the two constraints above. In 
the following, we briefly outline this procedure, focusing 
for simplicity on the half- filled case (n = 1/2) at zero- 
temperature. 

In this case, particle-hole symmetry makes the analysis 
simpler. It implies that (f> p = 4>n- p , and it can then be 
shown that the constraint X) p =i (^-i)^p = n = 1/2 is 
automatically satisfied when the first one (^2 p=0 (^)0p = 
1) is. Evaluating the ground-state energy of ([Ig|) yields 
the following expression to be minimized at T = (sub- 
ject to just one constraint): 



±E = 
N 



+ 



-4e 



'N-l 

E 

Lp=0 



N — 1 
P 



9 p+i 



TJ N 

-T 

2 ^ 

p=0 



N 
~2 



(50) 



with eo the (absolute value) of the non-interacting kinetic 
energy per orbital (as defined in (|36|)), and U = U/N, as 
defined above. 

It is clear from the Hamiltonian (^6|) that this slave- 
boson approximation leads to a picture of quasiparticles 
with no residual interactions. The quasiparticle weight 
reads: 



Z = z z = 




(51) 



while the chemical potential shift A in the previous sec- 
tion corresponds in the present context to the Lagrange 
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multiplier associated with the constraint (44a) (it van- 
ishes at half- filling) . 

We also note that the non-interacting limit is correctly 
reproduced, thanks to the square-root normalization fac- 
tors [^7j introduced in (ff3|). In that limit, Z = 1 and 
E = -Ne . 

In order to analyze the Mott transition, we note that 
within this approach, an insulating solution can be found 
for arbitrary strength of the coupling U. It corresponds 
to the trivial solution: 



<Pn/2 



N 
N/2 



K = 0(P^^) (52) 



which yields the T = occupancies in the atomic limit. 
This solution is always an extremum of the variational 
energy (corresponding to a boundary), but not always 
a minimum. Indeed, a minimum with a lower energy 
exists for U < U C 2- To find this critical coupling, we 
perturb around the trivial solution. Close to {7 c2 , only 
the occupancies of the states with charge N/2 ± 1 mat- 
ter, and we set: (j> 2 N/2 = (^ 2 ) + • • • , (f> N / 2 ±i = 6(j>, 



and <p p — ■ ■ ■ otherwise (where (• • • ) denote terms that 
can be neglected) . The energy difference with the trivial 
insulating solution reads: 

((jv/2) + Gv/2-1)) „ , . . 

AE = -4e V /Af / J —8<t> 2 (53) 



U ( ( N 



N 



2 V \N/2 - 1 / \N/2 + 1 



Hence, the critical coupling at which the metallic solution 
disappears reads: 



U C 2 = 16e - 



/JV-1\ S 

Kn/2) , N + 2 



( N W N \ 
VAT/2/ \N/2-l) 



= 4e n 



N 



(54) 



This value coincides with that obtained by Lu |24j us- 
ing the Gutzwiller approximation. Remarkably, it agrees 
with our exact result in the N — > 00 limit. Together 
with the considerations above, this suggests that such 
a slave boson approach becomes asymptotically exact, 
at low-energy and close to the critical coupling U C 2, in 
this limit. This estimate of the critical coupling is al- 
ready quite accurate for the smallest value of N — 2 (we 
take for reference a semi-circular band with half-width 
D = l): 

[/ C (/V = 2) = — (/V + 2)~3.3 (55) 

(the NRG and the projective analysis of the one-orbital 
model indicate that U c is very close to 3). On the other 
hand, the large- TV solution given by equation ( p5| ) misses 
obviously the sub- leading term of order TV , and is there- 
fore inaccurate at small N. 



The behaviour of the quasi-particle weight close to the 
transition can be obtained by expanding the energy to 
next order in 8(f), taking into account that 4^/ 2 ± q ~ (1 — 

U /U C 2) q (so that only the terms with q = 0,±1 matter). 
This yields: 



Z = 1 - 



u72 



(56) 



We note that there are no 1/N prefactors in this expres- 
sion when couplings are scaled properly. 

It is not entirely obvious to decide what the critical 
coupling U c i is within this approach. Since the insulat- 
ing solution exists for arbitrary U , one might be tempted 
to conclude that U& = 0. Another criterion would be 
to find at which coupling the optical gap of the insu- 
lator vanishes. Unfortunately, the above approximation 
(condensed bosons) is insufficient to discuss high-energy 
features (Hubbard bands). Studies that go beyond this 
approximation jp9[ and incorporate fluctuations of the 
slave bosons suggest that the optical gap closes at [7 c2 : in 
that sense the slave boson approximation does not yield 
a coexistence regime, in contrast to the exact results es- 
tablished above (which are the generic situation within 
DMFT). 



V. FINITE TEMPERATURE TRANSITIONS 

We conclude this paper by addressing the finite- 
temperature aspects of the transitions discussed above. 
As in the one-band case, we expect a coexistence re- 
gion U c i(T) < U < U C 2{T) at finite temperature, 
which closes at a second-order (liquid-gas) critical point 

(C/mit, Tmit), with J7 c i(^mit) = U C 2(Tmyt) = C^mit- 
Comparison of frcc-cncrgics then yield a first-order tran- 
sition line T = T C {U) for T < T MIT . This has been the 
subject of many studies in the single orbital case, and is 
now established on firm grounds. However, no analytical 
estimate of the critical temperature Tmit is available. 
This is highly desirable, since this scale appears to be 
strongly reduced in comparison to the electronic band- 
width. It would also be interesting to know how this 
scale depends on the orbital degeneracy. 

Let us start with a qualitative argument, which has 
two merits in our view: (i) it allows to understand why 
?mit A i s such a small energy scale for a small num- 
ber of orbitals and (ii) it provides a lower bound on 
?mit- This argument is based on a comparison of the 
energy difference between the metallic and insulating so- 
lutions at T = 0, to the corresponding entropy differ- 
ence. For U = NU close to U C 2, the energy difference 
AE = E r (T = 0) - E M (T = 0) is expected to behave 
as: 



AE = N 



(Uc2- Uf 
Uc2 



(57) 
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This is indeed supported by the slave-boson calculation 
of Sec. IV. The entropy difference AS = Si — Sm is dom- 
inated by the entropy of the insulating solution (which 
has a degenerate ground-state) since the entropy of the 
metal vanishes linearly with T at low temperature T. 
Therefore: 



AS ~ In 



U/2 



N In 2 



(58) 



This allows to estimate the behaviour of the first-order 
critical line U C (T) (or T C (U)) at low-temperature (i.e for 
U close to U C 2). Indeed, the transition into the insulating 
state upon raising temperature occurs when the entropic 
gain overcomes the energy cost, leading to: 



Tc(U) 



(U c2 - U f 



U, 



(59) 



c2 



This estimate is valid only at very low temperature, close 
to the T — critical point U — U C 2- It indicates a 
quadratic dependence of the critical line, with a prefactor 
which is independent of TV 30 1. 

When U c i is close to U C 2, which is the case for a 
small number of orbitals (in particular in the one-band 
case), one can roughly estimate the critical endpoint 
Tmit by using ( f59| ) in the regime of couplings where 
the insulating solution disappears, leading to: Tmit — 
(U C 2 — U c i/N) 2 /U C 2- It is clear from this formula that 
the ratio Tmit A i s small because of the small energy 
difference between the metallic and insulating solutions. 
For a larger number of orbitals however, this can only 
produce a lower bound on the critical endpoint Tmit- 
Indeed, the energy difference close to Tmit is certainly 
underestimated by (p7[), while the entropy difference is 
reduced as compared to A In 2 since the entropy of the 
metal cannot be neglected at Tmit (we note that it be- 
haves as 7 T with 7 oc N) . This argument shows that 
Tmit A mu st either saturate to a constant, or increase 
with N, as N increases. 

We now turn to a more quantitative study of these 
finite-temperature issues, using the Landau functional 
formalism introduced in JT4J] . A functional of the effec- 
tive hybridisation is introduced, which reads, for the N- 
orbital model on the Bethe lattice (semi-circular d.o.s): 



In Zi mp [A] 



(60) 



in which Zi mp is the partition function Z — Tr e imp of 
the effective impurity model, defined by the action: 



r0 

Simp = / dr d+d T d,j 
Jo 



U 
~2 



1>X 



f dr f dr'A(r - r') £ d+(rK(0 (61) 
Jo Jo 



This functional is locally stationary for these hybridisa- 
tion functions which satisfy the DMFT equations: 

SF 

— = «■ A(ito n ) = t 2 G{iuj n ) (62) 

The local stability of these DMFT solutions is controlled 
by the matrix of second derivatives: 



S 2 F 



SA(iuj n )5A(iuj r , 



N 

~2 ~l~ R-nm 



(63) 



with, using (pOjfirj): 



K nm = ^ /(f • • • Jo dndfidndri e -4-i-n)+*, m (^-r 2 ) 
£ CT1CT2 [( T r < (ri)d+ {r 2 )dM)) - 
(T r d+ ( n )d ai (r{)> (T T d+ (r 2 )d a2 (t' 2 ))\ (64) 

A solution of the DMFT equations is locally stable pro- 
vided the Xnm matrix has no negative eigenvalues when 
evaluated for this solution. Hence, the couplings U C \(T) 
(resp. U C 2{T)) correspond to the instability line of the 
(U,T) plane along which a negative eigenvalue appears 
when Xnm is evaluated for the insulating (resp. metal- 
lic) solution. For T = 0, these instability criteria should 
coincide with those derived from the projective method 
described above (as verified below for U C \(T = 0)). The 
critical endpoint T = Tmit is such that, for T > Tmit, 
no negative eigenvalue is found, at any coupling U , when 
the stability matrix is evaluated on the (unique) DMFT 
solution. It should thus be noted that, in order to de- 
termine ({/mit, Tmit), it is not necessary to know both 
U c i(T) and U C 2(T): either one of them is in principle 
sufficient. 

The practical difficulty in performing this stability 
analysis at finite temperature is that: (i) it requires a 
knowledge of the finite-T solution of the DMFT equa- 
tions -or at least a reasonable app roximation to it- and 
(ii) the two-particle correlator (^4() must be evaluated for 
this solution. Completing this program, for arbitrary or- 
bital degeneracy, using numerical methods, is beyond the 
scope of this paper. However, we would like to present 
here a simpler calculation which gives some insight in the 
dependence of f/ c i(T) (and of Tmit) on the orbital degen- 
eracy TV. What we have done is to use the atomic limit 
as a very rough first approximation to the insulating so- 
lution. We have evaluated the 2-particle correlator K nm 
and the stability matrix Xnm in this limit, for arbitrary 
N. This can be directly implemented on the computer, 
using the spectral decomposition of K nm onto atomic 
eigenstates (some details are provided in Appendix ^) . 
The stability matrix is then diagonalised numerically (a 
truncation to a -large- number of Matsubara frequencies 
is made), and we search for the line in the (U,T) plane 
where a negative eigenvalue is first found. 

The result of this calculation (for the half-filled case 
n = 1/2) is depicted in Fig. |l|. There, the temperature 
below which a negative eigenvalue is found is plotted as 
a function of U. We first observe that no instability is 
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found above a critical value of U 
with the value 



which does coincide 
4\/TV~ 



fl!6aD : U cl (T = 0)/t ~ 4\/TV + 1 esti- 
mated above from the projective method (in the atomic 
approximation). This is a non-trivial consistency check 
on our calculations. Secondly, we find that the insta- 
bility lines for different values of the orbital degeneracy 
TV can all be collapsed on a single curve when both T 
and U are rescaled by ^JN + 1. This can be understood 
from the fact that, in the atomic limit, only the two 
scales U and T enter the two-particle correlator K nm 
(not t), and that K is of order TV 2 , so that the instability 
criterion can be written (from dimensionality considera- 



tions): f = $(- 



: ). Naturally, we do not expect the 



curve in Fig. [l] to be a quantitatively reliable determina- 
tion of U c i(T), because of the atomic approximation in- 
volved. In particular, we see that, instead of terminating 
at a critical endpoint (J/mit, Tmit), the instability line in 
Fig. |l| bends back, yielding a stability window at small 
U which is certainly a spurious aspect of the approxima- 
tion. Hence only the branch of the curve corresponding 
to larger values of U has physical significance. 

It is tempting to conclude, from the the observed scal- 
ing with \/N + 1, that Tmit itself will grow in that man- 
ner as TV is increased. Indeed, we note that recent QMC 
calculations by Amadon and Biermann |2^| do show a 
marked increase of Tmit with TV, not inconsistent with 
a VTV scaling. However, we do not consider this issue 
to be entirely settled: rather, this atomic estimate pro- 
vides an upper bound on the growth of Tmit with TV. 
Indeed, it neglects the presence of the quasiparticle reso- 
nance (which is reduced but still present as temperature 
is raised close to Tmit)- The resonance makes the solu- 
tion less stable, so that the atomic calculation is likely to 
overestimate Tviit- Combined with the above entropy ar- 
gument (which provides a lower bound), we conclude that 
Tmit A i s an increasing function of TV at moderate val- 
ues of TV, which either saturates at large- TV or increases 
at most as y/~N. Clearly, a more refined analysis of the 
above stability condition is required to settle this issue. 



VI. CONCLUSION 

In this paper, we have studied the Mott transition 
of the TV-orbital Hubbard model in the fully symmetric 
case. The physical picture is qualitatively the same as for 
TV = 2, i.e. a coexistence region exists between an insu- 
lating and a metallic solution within two distincts critical 
couplings Uci and U c %- We have shown that these critical 
couplings do not have the same dependence on the num- 
ber of orbitals, as the latter becomes large: U c \ cx \/~N 
while U C 2 cx TV. We have obtained an exact analytical 
determination of the critical value U C 2 — U c ijN in this 
limit. These results explain the widening of the coexis- 
tence window \U C \ , U C 2\ observed in the numerical simula- 
tions of Ref . p5] , and more recently in p^] . Our findings 
also put the results of Refs. [H], in a new perspective, 
as discussed in the introduction. 



T 0.04 




FIG. 1: Instability temperature below which a negative eigen- 
value is found when the stability matrix Xnm of Eq. (p3|) is 
evaluated in the atomic approximation. The various curves 
are for increasing number of orbitals (TV = 2,4,6,8,10,20), 
and an almost perfect scaling with T/vTV+1 and U/^/N + 1 
is found. 



We have shown that, at large TV and close to the crit- 
ical coupling U C 2, the scaling U — NU is appropriate, so 
that the separation of scale which is at the heart of the 
projective method becomes asymptotically exact. The 
low-energy physics close to U c i is then well described by 
slave-boson like approximations. This is the case in par- 
ticular for the quasi-particle peak in the spectral func- 
tion, which is separated from Hubbard bands by a large 
energy (of order NU) and is accurately described by the 
simple form (|32|). In view of practical applications to 
electronic structure calculations with large orbital degen- 
eracy (e.g f-electron systems), it would be highly desir- 
able to have also an analytical determination of the high- 
energy part of the spectral function in the large- TV limit 
(or at least reasonable approximations to it). This is left 
for future investigations. 

Another issue that we have only partly addressed is 
the finite-temperature aspects of these transitions. The 
whole finite-temperature coexistence region widens as TV 
is increased, and the critical temperature associated with 
the Mott critical endpoint does increase with the num- 
ber of orbitals at intermediate values of TV. The precise 
behaviour of Tmit A at large TV deserves further studies 
however. 
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APPENDIX A: LARGE-TV SCALINGS AND THE 
ATOMIC LIMIT 

In this appendix, we provide a detailed proof that the 
corrections in Vh to the atomic limit (Vh = 0) do not 
modify the large- N dependence of the critical couplings. 
This is done by an explicit investigation of the struc- 
ture of these corrections. Moreover, we show that for 
U oc N (i.e when U C 2 is considered), these corrections 
only produce subdominant terms which can be ignored 
at large- N. 



Corrections in Vh for U = NU. 



For our purpose, the operator ifjj can be written : 



K 



Hi 



K 



(Al) 



U 
2~ 



E 



feeH±,cr 



{d) a d a - n) 



First let us consider the ground state energy E gs . As 
the lower Hubbard band H - is filled in \n) , we prefer to 
remove the constant contribution X)fceH- a so that the 
zeroth order energy is Eo = o(a*| H^ Im)o = 0- The next 
order term <5Ei = o(mI Ia*)o * s zer0 by conservation of 
the total charge Q on the d-level, so we need to look at 
the second order contribution: 



<5E 2 = £ 



Eo — E 



(A2) 



= -(N-Q) 



V 2 
1 H 



fceH 



- 2 ~ £ fe 



1 E 

fe6H+ 



V 2 



T + £ fe 



since |V>) can only possess a charge Q + 1 (resp. Q — 1) 
for the d-electron and a hole (resp. particule) excita- 
tion in the bath. Clearly as tk is negative (positive) for 
k € H~ (H + ), the denominators are of order N (since 
U = NU) and <5E 2 is therefore scaling as N° at large 
N. We would like to infer this to be true at all order in 
the developement in Vh- One can easily classify into two 
categories the intermediate states that appear at order 
V[f and beyond. There can either be a charge different 
form Q on the impurity (as in the previous computation 
at order ) and then the denominator coming from this 
state provides (at least ) a factor TV -1 by the presence 
of the Coulomb energy U = NU. Or the impurity can 
be found in a charge Q state (with possibly a reorganiza- 
tion of the spin configuration) without a Coulomb energy 
cost. However this state possesses also many particle- 
hole excitations in the bath, and this leads to a denom- 
inator which scales as the Mott gap A g . For example 
let us consider such a ket generated at order V£, like 



| <7i . . . <tq) d |2 . . . N)_ 1 1) + . This particular state (but the 
argument is general), contributes a factor (e^ — £fc<) _1 
where k G H~ and k' € H + . This quantity is always 
smaller than A" 1 because the Mott gap A g is also the 
gap of the Ck,v electrons (through the self-consistency). 
As the gap A g is expected to be of order U = NU when 
U is close to U C 2, each particle-hole excitation also pro- 
vides a factor iV _1 in the perturbative expansion. This is 
enough to conclude that, at a given order in perturbation 
theory, each of the denominator gives a factor TV -1 that 
balance the combinatorial factor coming from the string 
of charge excitation on the impurity level. Thus E gs is 
at most of order 1 when N is large. 

We finally need to compute J™ at large N. Let us 
consider first the case fj, ^ v, so that /i and v will differ 
by only one spin flip a <-» a'. We can now calculate at 
the lowest order in Vn the desired matrix element: 



(n\d a 



1 



K 



Hi 



E n 



(A3) 



with 



Im) = Im) + Im)i 



(A4) 



Hi = E v H 4 c fc,gl/*)o + E u Vh . c L^Im) 



u 

fcGH-,o- 2 



keH+,a 2 + ( 'k 



The contribution at order (Vh)° is therefore 



'B^*M*~V (A5) 



because E gs = 0(1) and U = NU. All terms coming 
at order Vh cancel by conservation of the charge, so we 
examine now the next leading order for T, ST 2, which 
is composed of three terms: one from the correction by 
Vh to each of the two external k ets, one from the de- 
velopment of the denominator in (A5) at order V£ , and 
a mixed term (between only one ket correction and the 
denominator developed at first order). Let us examine 
the first of these second order contributions: 



ST 



(i) 



1 



H u - E gs 



dtWh 



(A6) 



The state is composed of charge Q + 2 excitations 

(with energy 2 2 U/2) and of particle- hole excitations in 
the bath with a charge Q on the impurity (this has the 
energy Ck — £*/)• Both energies scale like N and con- 
tribute an overall factor N~ 3 to ST^ ■ There is also a 
combinatorial factor of order N coming from the choice 
in the spin-excitation in but there is none from the 
state I /j,) 1 as the spin configuration in d a , |/i)i is fixed by 
the one in d\, \v) x . Therefore ST^ = 0(N- 2 ). The same 
argument works for the two other contribution at order 
V 2 

More generally, at a given order p in the development 
in Vh, each of the p+l denominators provides a factor of 
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order N 1 , because there is a finite gap A g ~ f/ C 2 cx N this development provides a leading contribution of order 



between the two Hubbard bands 1 25 . However the string 
of p spin-flips generated by the H mlx terms only con- 
tributes to an overall combinatorial factor at most N p / 2 
because one has to connect two fixed external configura- 
tions \v) and Hence in the limit of large N all terms 
after the first one are sub-dominant, so that: 



r = <5r n 



o 



i 



2 
U 



2 
NU 



(A7) 



meaning that the atomic result for T is exact in the large- 
TV limit when U cx N. 

At this stage, we want to stress a posteriori that the 
hypotheses made in the preceding argument are consis- 
tent with this result. First it is correct to fix the chemical 
potential Jl — from the atomic limit (Vh = 0), because 
perturbative corrections do not change the occupancy of 
the d-orbital, as can be checked from a direct expansion 
in Vh of (52 a d^dcr). Finally the existence of a large Mott 
gap A g cx N can also be justified on the same grounds. 



1/y/N ( using U ~ A 9 cx VN). T can not be computed 



exactly in that case, because B(N, n; Vn/(V NU)) is now 
a non trivial function of Vr/U at N = oo. But this is 
enough to conclude that B(N,n;Vn/U) is of order 1 at 
large N when U cx \/N, so that the result U c \ ~ \/N~U c \ 
holds. However, an exact evaluation of U c \ appears to be 
a difficult task, even at N = oo. 



APPENDIX B: NUMERICAL CALCULATION OF 
THE STABILITY MATRIX IN THE ATOMIC 
LIMIT 



In this Appendix, we provide details on the computa- 
tion of the stability matrix Xnm in the A-orbital atomic 
limit (cf. Fig. [J). First, we use a spectral decomposition 
over the atomic eigenstates in Eq. (f34|). This yields: 

K nm = A(iu n , ~iuj n , iv n , -iv n ) / Z at - B(ioj n )B(iv n )/Z% t 



Some remarks about U c i 



When U cx yJV on the other hand, the perturbative 
series for T does not appear to stop after the first contri- 
bution, showing that U c \ is affected by Vh in the large A" 
limit. It is however straigh tforw ard to follow the line of 
arguments given in Section [II C to see that each term in 



with 



N 

Zat = ]T e"^« 

Q=0 



Afa,j = 1..4) e"^« e (P)(d(Q,P) + (AT- l)C 2 (g,P))p 4 (wp W +A Q+Ef= . +iaj5(i))0p( . ) ,j = 1..4) 



0<Q<N 
P£Sa 



N 



B{iu n ) =NJ2 



-0Eq 



Q=0 



N -I 
Q-l 



F 2 (iu n + Aq-i,!, -iu n + Aq : _i) - 



A- 1 

Q 



F 2 (-iu: n + AQ+^-ijiWn + A Qi i) 



'< P ){q-2) if P" 1 (l)<^ 1 (2) and P~\3) < P" x (4) 
C 2 (Q,P)= \e(P)( N Q 2 ) if P-1(1) >p-\2) and i>-*(3) > P- x (4) 
-e(P)( N Q Zl) else 



Ci(Q,P) 



< ,N Zl) if P- 1 ({1,3}) = {1,3} 
'V) if P- 1 ({1,3}) = {2,4} 



else 



F n (uj,j = l..n) 



/3>Ti>->T„ 



»=i 



In these expressions, Eq = U(Q — N/2) 2 /2 are the 
atomic energy levels, Aq )0 = Pg+ a — Eq, e is the sig- 
nature of the permutation P, and we used the nota- 
tion f(xj,j = 1..4) = f(xx, x 2 , X3, X4). We compute P4 
and F2, using the relations io\ + u>2 + W3 + uj± = and 



u)\ + u>2 = respectively. The algorithm can be decom- 
posed into three functions : (i) From U and T, compute 
Xnm for |n|, |m| < 100 and then its lowest eigenvalue 
Eo(U,T); (ii) Find a point pO below the instability tem- 
perature Ti and a majoration Tm of max(Ti) and [7^ 
of C/ C i; define a path in the (U,T) from (0,0) to (0,T M ) 
to (Um,Tm) to (£7jw,0); (iii) For about 50 points p on 
this path, solve for Eq(U, T) — on the line defined by p 
and po using a dichotomy. To speed up the computation, 
the sums over permutations and Q in A and B are ex- 
panded automatically into C++ code, which is inlined in 
the main program. Diagonalisations are performed using 
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the Fortran LAPACK library. Codes will be available at www-spht.cea.fr/~parcolle/. 
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